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ABSTRACT 

We examine several recent calculations of the self-energy of 
a neutron star arising from neutrino-exchange. It is shown that 
the results of Abada, et al. in 1+1 dimensions have no bearing 
on a 3-dimensional neutron star, since the criticality parameter 
GfN/R 2 is always much smaller than unity in 1+1 dimensions. 
The calculation of Kiers and Tytgat in 3-dimensions is shown 
to disagree with the lowest order 2-body contribution, which is 
known exactly. This discrepancy raises the possibility that the 
description of a neutron star as a continuous medium may be 
inappropriate when calculating higher-order many-body effects. 
We conclude that none of the recent calculations contradict the 
earlier claims that the neutrino-exchange contributions to the self- 
energy of a neutron star are unphysically large, when calculated in 
the standard model. The implication of this result, that neutrinos 
must have a non-zero mass, m u > 0.4 eV/c 2 , remains intact. 



It has been argued that many-body forces arising from the exchange of 
massless neutrinos can lead to an unphysically large energy density in white 
dwarfs and neutron stars [1, 2]. Intuitively this comes about because in the 
absence of any other physical parameters, the /c-body (k = 2, 4, 6, ...) 
contribution to the binding energy W of a neutron star of radius R 

must be of the form [1, 2] 
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Here G F is the Fermi constant and (^j is the binomial coefficient, 
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Combining Eqs. (1) and (2) we find [1, 2] 



for N->k. (2) 
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For a typical neutron star (G F N/R 2 ) = O(10 13 ), from which it follows that 
for k <C N higher-order many-body interactions make increasingly larger 
contributions to [3]. The energy, W = Yi k W^ can exceed the mass of 
the neutron star, as is shown explicitly in Ref . [1] , and this eventually leads to 
the conclusion that neutrinos must have a minimum mass, m„ > 0.4 eV/c 2 . 
The effect of a nonzero m v is to produce a "saturation" of the neutrino- 
exchange force, and ultimately a physically acceptable value for W. 

Following the publication of Ref. [1] , a number of papers appeared dealing 
with various questions raised by this calculation [4-10]. Here we focus on 
attempts to calculate the energy density w non-perturbatively starting from 
variants of the Schwinger formula [11] used in Ref. [1]: 
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In Eq. (4) a n = — 1/2 is the neutrino-neutron coupling constant, is the 
neutrino current, and Sp\E) is the free neutrino propagator. In the no- 
tation of Ref. [1], = -ip(x)(5 M 4, where p[x) is the neutron density, with 
/ d 3 xp(x) = N = O(10 57 ). The perturbative method of Ref. [1] treats neu- 
trons as discrete particles, and W can be calculated for any p{x). Our results, 
to be presented elsewhere show that the result in Ref. [1] is robust: For a va- 
riety of symmetric and asymmetric matter distributions the calculation of W 
given in Ref. [1] leads to an unphysically large value for massless neutrinos. 

By constrast the non-perturbative calculations treat the neutron star as 
a continuous medium whose density is constant, either throughout all space 
[4, 8] or in some sub-region [7, 10]. Since various cancellations occur in both 
approaches, one explanation for the differences among the results may be 
the different pictures used to describe a neutron star (i.e. discrete neutrons 
versus a continuous medium), as we discuss in more detail elsewhere. Along 
with other considerations, treating the neutron star as an object of finite 
extent is essential. 

Two recent papers have attempted to calculate W for a "neutron star" of 
finite extent, rather than calculating the energy density of an infinite medium. 
Abada, et al. [7] have calculated the effect of many-body neutrino exchange 
in a 1+1 dimensional "neutron star", and Kiers and Tytgat [10] have carried 
out a similar calculation in 3+1 dimensions. 

Returning to Eqs. (l)-(3) we note from the heuristic derivation of Eq. (3) 
- or from the formal treatment in Ref. [1] - that will have the same 
functional form in any number of spatial dimensions. The differences among 
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the 1-, 2-, and 3-dimensional cases arise when (GfN/R 2 ) is re-expressed in 
terms of the corresponding linear, surface and volume densities: 

1-dimension : For a linear array of length L and density A, N = XL 
and R = L/2 so that 

G F N/R 2 = AG F X/L . (5) 

The maximum value of A corresponds to taking the inter-neutron spacing 
to be given by the hard core radius r c = 0.5 x 10 -13 cm, which gives A = 
2 x 10 13 cm" 1 . Taking R = R 1Q = 10 km then gives 

G F N/R 2 10 \ 1 _ dim = 2x 10- 25 . (6) 

Hence in 1-dimension the "criticality parameter" G F N/R 2 in Eq. (3) is 
always much less than unity, from which one could conclude at the outset 
that W would always remain acceptably small. Physically this comes about 
because in 1-dimension there are too few neutrons in the neighborhood of a 
given neutron to allow the combinatoric factor (J^J to provide a sufficiently 
large enhancement. 

From Eq. (6) we recognize that many-body (k > 2) effects are supressed 
in 1-dimension, and hence the leading contribution to the energy comes from 
the 2-body neutron-neutron potential [1, 2, 12-16], 

r 2 n 2 
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The average interaction energy U{_ dim of a pair of neutrons in a linear array 
of length L would be given by 
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Where Pi(r) is the 1-dimensional probability density given by [1] 

2 2r 

«(r> = §-£ (9) 

The leading term to the average energy is found by combining Eqs. (8) and 

(9), 

(2) _ ^ F a n , . 

We note that this result is unaffected by inclusion of any contributions from a 
trapped neutrino sea, since the 2-body potential is dominated by small values 
of r [9]. For a linear array of N discrete neutrons, there are approximately 
N 2 /2 pairs that can be formed, hence the energy in 1-dimension is 

= W<i> lm = 3fj££. („) 

Consequently, the energy density u>i_di m is a constant, given by 

— . = ^. (12) 

This result conflicts with the 1-dimensional results of Abada et al. [7] who 
find tui-dim to be identically zero. We note that since the Abada et al. 
results can be expanded in a power series in Gf, they must agree with the 
lowest-order perturbative results given in Eqs. (11) and (12). 

2-dimensions : For a planar circular array of neutrons with surface 
density a, we have iV = airR 2 so that 

G F N/R 2 = irG F v = constant. (13) 

Assuming, as before, a hard core radius r c = 0.5 x 1CT 13 cm gives a = 
1 x 10 26 cm -2 , and for R = R w = 10 km we find 
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G F N/R 2 w \ 2 _ dim = 2x l(T b . (14) 

Since this is again smaller than unity, we conclude that the neutrino-exchange 
energy W is well-behaved in 2-dimensions and dominated by the 2-body 
contribution. The energy density in 2-dimensions u>2-dim, can be obtained in 
a similar manner to the 1-dimensional result of Eq. (12), and we find 

r 2 n 2 rr 2 

«*-*- = 1W' (15) 
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3-dimensions : For a spherical volume of radius R and number 
density p we have in 3-dimensions iV = |7i\R 3 p, which gives 

G F N/R 2 = (4/3)7rG F pR. (16) 
Using p = 4 x 10 38 cm -3 from Ref. [1] gives 

G F N/R 2 10 \ 3 _ dim = 8 x 10 12 , (17) 

and hence it is only in 3 spatial dimensions that the possibility exists for 
an unphysically large energy density arising from the exchange of massless 
neutrinos. Thus the recent calculation of Abada, et ai, which is done in 
1-dimension, has little bearing on the calculation in Ref. [1]. 

We turn next to the paper of Kiers and Tytgat (KT) [10] which also 
calculates W under the assumption that a neutron star can be described as 
a continuous medium. KT work in a finite 3-dimensional volume, and hence 
their starting point is closest to that of Ref. [1] amongst the recent papers 
which use the continuous medium approximation. Since KT do not find a 
large value for W, in contrast to Ref. [1], it is instructive to compare the 
discrete and continuum calculations of W. 
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For present purposes we focus on the 2-body contribution to the 
self-energy of a neutron star W, which can be calculated exactly with no 
approximations [2]. The 2-body potential is given by Eq. (7). The aver- 
age interaction energy C/( 2 ) of a pair of neutrons having a uniform density 
distribution in a spherical volume of radius R is then given by [1, 2] 

where P(r) is the probability density function^ given by [1, 2] 

3r 2 



! 3 / r \ 1 / r \ 3 
~ 2 V2~rJ + 2 V2R~) 



v(r c ,R), (19a) 



i2) = --3 --e • (19b) 

l-8fe) +9fe) -2(&) 

Combining Eqs.(18) and (19) we find that the energy per pair of neutrons, 

U^ 2 ' is exactly given by 

^s^^-S) 3 " <*> 

where we have reinstated the factor he. For a neutron star containing N dis- 
crete neutrons, there are N(N — 1)/2 pairs and hence the 2-body contribution 
to W is given by 



(2) = 3 (G F q n ) 2 iV(iV-l) / _ r^N 
" 16tt 3 he B?rl V 2i?y 



%(r c ,i?), (21) 



which is also an exact result. Since iV 3> 1 for neutron star, we can replace 
N(N — 1) by N 2 , and re-express N in terms of the number density p: 
lr rhe normalization factor rj(r c ,R) was approximated by unity in Refs. [1] and [2]. 
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rj(r c ,R). 



(22) 



Also of interest is the energy density = /V, 



w 



1 (G F a n p) 2 
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r)(r c ,R). 



(23) 



It follows from Eq. (23) that the energy density is an increasing function 
of R, and approaches a constant as R — > oo. 

We now compare these results to those obtained by KT for the 2-body 
contribution. KT find for 



where A is a cutoff which plays a role similar to r c in Eq. (21). KT then use 
Eq. (24) to conclude that 



Comparing Eqs.(24) and (25) to Eqs.(22) and (23), respectively, we see that 
the independence of the KT results disagrees with that of the exact results 
obtained by assuming discrete neutrons. Since the KT results can be ex- 
panded in a power series in Gp, the lowest order KT results must agree 
with the corresponding perturbative results given in Eqs. (22) and (23). 
Furthermore the KT results are seen to underestimate and the energy 
density as R increases. This is a particularly troublesome point since 
it raises the possibility that the many-body contributions may be similarly 
underestimated in the KT formalism. 




(24) 




(25) 
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It should be emphasized that the discrepancy between the KT results 
and the exact results is a matter of principle, and cannot be dismissed on the 
grounds that is in either case small. Until the origin of this discrepancy 
is fully understood one cannot be certain that it does not affect other parts 
of their calculation as well. 

CONCLUSIONS 

The recent papers by Abada, et al. [7], and by Kiers and Tytgat [10] 
have been analyzed. We have shown explicitly that the results of Ref. [7] 
in 1-dimension have no bearing on whether the neutrino-exchange energy 
density in 3-dimensions is unphysically large in a neutron star. The KT cal- 
culation, although carried out in 3-dimensions, fails to reproduce the 2-body 
result which is known exactly and, moreover, underestimates this contribu- 
tion. Obtaining the correct 2-body result is a litmus test for the validity of 
any calculational scheme, and hence the significance of the KT calculation 
is unclear at present. One question that must be explored is the validity 
of approximating the neutron star by a continuous medium. Although such 
a picture is routinely employed in discussing the MSW effect (in which Gf 
enters only in lowest order), its appropriateness for calculating self-energy in- 
teractions involving higher-order weak processes remains to be demonstrated. 
Another question that remains to be understood is where the medium de- 
scription breaks down: Clearly an atomic nucleus should be described in 
terms of discrete neutrons for purposes of calculating W . Since a neutron 
star behaves in some way as a big nucleus, the transition from a description 
in terms of discrete neutrons to one in terms of a medium must be fully un- 
derstood. These and other related questions will be explored in more detail 
elsewhere. 
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